GWAS reveals determinants of mobilization rate and dynamics of an active endogenous retrovirus of cattle.

Five to ten percent of mammalian genomes is occupied by multiple clades of endogenous retroviruses (ERVs), that may count thousands of members. New ERV clades arise by retroviral infection of the germline followed by expansion by reinfection and/or retrotransposition. ERV mobilization is a source of deleterious variation, driving the emergence of ERV silencing mechanisms, leaving "DNA fossils". Here we show that the ERVK[2-1-LTR] clade is still active in the bovine and a source of disease-causing alleles. We develop a method to measure the rate of ERVK[2-1-LTR] mobilization, finding an average of 1 per ~150 sperm cells, with >10-fold difference between animals. We perform a genome-wide association study and identify eight loci affecting ERVK[2-1-LTR] mobilization. We provide evidence that polymorphic ERVK[2-1-LTR] elements in four of these loci cause the association. We generate a catalogue of full length ERVK[2-1-LTR] elements, and show that it comprises 15% of C-type autonomous elements, and 85% of D-type non-autonomous elements lacking functional genes. We show that >25% of the variance of mobilization rate is determined by the number of C-type elements, yet that de novo insertions are dominated by D-type elements. We propose that D-type elements act as parasite-of-parasite gene drives that may contribute to the observed demise of ERV elements.

Sequenced-based GWAS for linear classification traits in Belgian Blue beef cattle reveals new coding variants in genes regulating body size in mammals.

BACKGROUND: Cohorts of individuals that have been genotyped and phenotyped for genomic selection programs offer the opportunity to better understand genetic variation associated with complex traits. Here, we performed an association study for traits related to body size and muscular development in intensively selected beef cattle. We leveraged multiple trait information to refine and interpret the significant associations. RESULTS: After a multiple-step genotype imputation to the sequence-level for 14,762 Belgian Blue beef (BBB) cows, we performed a genome-wide association study (GWAS) for 11 traits related to muscular development and body size. The 37 identified genome-wide significant quantitative trait loci (QTL) could be condensed in 11 unique QTL regions based on their position. Evidence for pleiotropic effects was found in most of these regions (e.g., correlated association signals, overlap between credible sets (CS) of candidate variants). Thus, we applied a multiple-trait approach to combine information from different traits to refine the CS. In several QTL regions, we identified strong candidate genes known to be related to growth and height in other species such as LCORL-NCAPG or CCND2. For some of these genes, relevant candidate variants were identified in the CS, including three new missense variants in EZH2, PAPPA2 and ADAM12, possibly two additional coding variants in LCORL, and candidate regulatory variants linked to CCND2 and ARMC12. Strikingly, four other QTL regions associated with dimension or muscular development traits were related to five (recessive) deleterious coding variants previously identified. CONCLUSIONS: Our study further supports that a set of common genes controls body size across mammalian species. In particular, we added new genes to the list of those associated with height in both humans and cattle. We also identified new strong candidate causal variants in some of these genes, strengthening the evidence of their causality. Several breed-specific recessive deleterious variants were identified in our QTL regions, probably as a result of the extreme selection for muscular development in BBB cattle.

The rate of de novo structural variation is increased in in vitro-produced offspring and preferentially affects the paternal genome.

Assisted reproductive technologies (ARTs), including in vitro maturation and fertilization (IVF), are increasingly used in human and animal reproduction. Whether these technologies directly affect the rate of de novo mutation (DNM), and to what extent, has been a matter of debate. Here we take advantage of domestic cattle, characterized by complex pedigrees that are ideally suited to detect DNMs and by the systematic use of ART, to study the rate of de novo structural variation (dnSV) in this species and how it is impacted by IVF. By exploiting features of associated de novo point mutations (dnPMs) and dnSVs in clustered DNMs, we provide strong evidence that (1) IVF increases the rate of dnSV approximately fivefold, and (2) the corresponding mutations occur during the very early stages of embryonic development (one- and two-cell stage), yet primarily affect the paternal genome.

An organism-wide ATAC-seq peak catalog for the bovine and its use to identify regulatory variants.

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

High-resolution structural variants catalogue in a large-scale whole genome sequenced bovine family cohort data.

BACKGROUND: Structural variants (SVs) are chromosomal segments that differ between genomes, such as deletions, duplications, insertions, inversions and translocations. The genomics revolution enabled the discovery of sub-microscopic SVs via array and whole-genome sequencing (WGS) data, paving the way to unravel the functional impact of SVs. Recent human expression QTL mapping studies demonstrated that SVs play a disproportionally large role in altering gene expression, underlining the importance of including SVs in genetic analyses. Therefore, this study aimed to generate and explore a high-quality bovine SV catalogue exploiting a unique cattle family cohort data (total 266 samples, forming 127 trios). RESULTS: We curated 13,731 SVs segregating in the population, consisting of 12,201 deletions, 1,509 duplications, and 21 multi-allelic CNVs (> 50-bp). Of these, we validated a subset of copy number variants (CNVs) utilising a direct genotyping approach in an independent cohort, indicating that at least 62% of the CNVs are true variants, segregating in the population. Among gene-disrupting SVs, we prioritised two likely high impact duplications, encompassing ORM1 and POPDC3 genes, respectively. Liver expression QTL mapping results revealed that these duplications are likely causing altered gene expression, confirming the functional importance of SVs. Although most of the accurately genotyped CNVs are tagged by single nucleotide polymorphisms (SNPs) ascertained in WGS data, most CNVs were not captured by individual SNPs obtained from a 50K genotyping array. CONCLUSION: We generated a high-quality SV catalogue exploiting unique whole genome sequenced bovine family cohort data. Two high impact duplications upregulating the ORM1 and POPDC3 are putative candidates for postpartum feed intake and hoof health traits, thus warranting further investigation. Generally, CNVs were in low LD with SNPs on the 50K array. Hence, it remains crucial to incorporate CNVs via means other than tagging SNPs, such as investigation of tagging haplotypes, direct imputation of CNVs, or direct genotyping as done in the current study. The SV catalogue and the custom genotyping array generated in the current study will serve as valuable resources accelerating utilisation of full spectrum of genetic variants in bovine genomes.

Transcriptome analysis reveals tumor microenvironment changes in glioblastoma.

A better understanding of transcriptional evolution of IDH-wild-type glioblastoma may be crucial for treatment optimization. Here, we perform RNA sequencing (RNA-seq) (n = 322 test, n = 245 validation) on paired primary-recurrent glioblastoma resections of patients treated with the current standard of care. Transcriptional subtypes form an interconnected continuum in a two-dimensional space. Recurrent tumors show preferential mesenchymal progression. Over time, hallmark glioblastoma genes are not significantly altered. Instead, tumor purity decreases over time and is accompanied by co-increases in neuron and oligodendrocyte marker genes and, independently, tumor-associated macrophages. A decrease is observed in endothelial marker genes. These composition changes are confirmed by single-cell RNA-seq and immunohistochemistry. An extracellular matrix-associated gene set increases at recurrence and bulk, single-cell RNA, and immunohistochemistry indicate it is expressed mainly by pericytes. This signature is associated with significantly worse survival at recurrence. Our data demonstrate that glioblastomas evolve mainly by microenvironment (re-)organization rather than molecular evolution of tumor cells.

A High-Quality Genome Assembly of Striped Catfish (Pangasianodon hypophthalmus) Based on Highly Accurate Long-Read HiFi Sequencing Data.

The HiFi sequencing technology yields highly accurate long-read data with accuracies greater than 99.9% that can be used to improve results for complex applications such as genome assembly. Our study presents a high-quality chromosome-scale genome assembly of striped catfish (Pangasianodon hypophthalmus), a commercially important species cultured mainly in Vietnam, integrating HiFi reads and Hi-C data. A 788.4 Mb genome containing 381 scaffolds with an N50 length of 21.8 Mb has been obtained from HiFi reads. These scaffolds have been further ordered and clustered into 30 chromosome groups, ranging from 1.4 to 57.6 Mb, based on Hi-C data. The present updated assembly has a contig N50 of 14.7 Mb, representing a 245-fold and 4.2-fold improvement over the previous Illumina and Illumina-Nanopore-Hi-C based version, respectively. In addition, the proportion of repeat elements and BUSCO genes identified in our genome is remarkably higher than in the two previously released striped catfish genomes. These results highlight the power of using HiFi reads to assemble the highly repetitive regions and to improve the quality of genome assembly. The updated, high-quality genome assembled in this work will provide a valuable genomic resource for future population genetics, conservation biology and selective breeding studies of striped catfish.

A 2-month field cohort study of SARS-CoV-2 in saliva of BNT162b2 vaccinated nursing home workers.

Background: Nursing home (NH) residents have been severely affected during the COVID-19 pandemic because of their age and underlying comorbidities. Infection and outbreaks in NHs are most likely triggered by infected workers. Screening for asymptomatic NH workers can prevent risky contact and viral transmission to the residents. This study examined the effect of the BNT162b2 mRNA COVID­19 (Comirnaty®; BioNTech and Pfizer) vaccination on the saliva excretion of SARS-CoV-2 among NH workers, through weekly saliva RT-qPCR testing. Methods: A 2-month cohort study was conducted among 99 NHs in the Walloon region (Belgium), at the start of February 2021. Three groups of workers, i.e., non-vaccinated (n = 1618), one-dosed vaccinated (n = 1454), and two-dosed vaccinated (n = 2379) of BNT162b2 mRNA COVID­19 vaccine, were followed-up weekly. Their saliva samples were used to monitor the shedding of SARS-CoV-2. All positive samples were sequenced and genotyped to identify the circulating wild-type virus or variants of concern. Results: The protection fraction against the excretion of the SARS-CoV-2 in the saliva samples of the workers after the second dose is estimated at 0.90 (95% CI: 0.18; 0.99) at 1 week and 0.83 (95% CI: 0.54; 0.95) at 8 weeks. We observe more circulating SARS-CoV-2 and a greater variability of viral loads in the unvaccinated group compared to those of the vaccinated group. Conclusions: This field cohort study advances our knowledge of the efficacy of the mRNA BNT162b2 COVID-19 vaccine on the viral shedding in the saliva specimens of vaccinated NH workers, contributing to better decision-making in public health interventions and management.

ABO genotype alters the gut microbiota by regulating GalNAc levels in pigs.

The composition of the intestinal microbiome varies considerably between individuals and is correlated with health1. Understanding the extent to which, and how, host genetics contributes to this variation is essential yet has proved to be difficult, as few associations have been replicated, particularly in humans2. Here we study the effect of host genotype on the composition of the intestinal microbiota in a large mosaic pig population. We show that, under conditions of exacerbated genetic diversity and environmental uniformity, microbiota composition and the abundance of specific taxa are heritable. We map a quantitative trait locus affecting the abundance of Erysipelotrichaceae species and show that it is caused by a 2.3 kb deletion in the gene encoding N-acetyl-galactosaminyl-transferase that underpins the ABO blood group in humans. We show that this deletion is a ≥3.5-million-year-old trans-species polymorphism under balancing selection. We demonstrate that it decreases the concentrations of N-acetyl-galactosamine in the gut, and thereby reduces the abundance of Erysipelotrichaceae that can import and catabolize N-acetyl-galactosamine. Our results provide very strong evidence for an effect of the host genotype on the abundance of specific bacteria in the intestine combined with insights into the molecular mechanisms that underpin this association. Our data pave the way towards identifying the same effect in rural human populations.

Reliable and Scalable SARS-CoV-2 qPCR Testing at a High Sample Throughput: Lessons Learned from the Belgian Initiative.

We present our approach to rapidly establishing a standardized, multi-site, nation-wide COVID-19 screening program in Belgium. Under auspices of a federal government Task Force responsible for upscaling the country's testing capacity, we were able to set up a national testing initiative with readily available resources, putting in place a robust, validated, high-throughput, and decentralized qPCR molecular testing platform with embedded proficiency testing. We demonstrate how during an acute scarcity of equipment, kits, reagents, personnel, protective equipment, and sterile plastic supplies, we introduced an approach to rapidly build a reliable, validated, high-volume, high-confidence workflow based on heterogeneous instrumentation and diverse assays, assay components, and protocols. The workflow was set up with continuous quality control monitoring, tied together through a clinical-grade information management platform for automated data analysis, real-time result reporting across different participating sites, qc monitoring, and making result data available to the requesting physician and the patient. In this overview, we address challenges in optimizing high-throughput cross-laboratory workflows with minimal manual intervention through software, instrument and assay validation and standardization, and a process for harmonized result reporting and nation-level infection statistics monitoring across the disparate testing methodologies and workflows, necessitated by a rapid scale-up as a response to the pandemic.

Benchmarking phasing software with a whole-genome sequenced cattle pedigree.

BACKGROUND: Accurate haplotype reconstruction is required in many applications in quantitative and population genomics. Different phasing methods are available but their accuracy must be evaluated for samples with different properties (population structure, marker density, etc.). We herein took advantage of whole-genome sequence data available for a Holstein cattle pedigree containing 264 individuals, including 98 trios, to evaluate several population-based phasing methods. This data represents a typical example of a livestock population, with low effective population size, high levels of relatedness and long-range linkage disequilibrium. RESULTS: After stringent filtering of our sequence data, we evaluated several population-based phasing programs including one or more versions of AlphaPhase, ShapeIT, Beagle, Eagle and FImpute. To that end we used 98 individuals having both parents sequenced for validation. Their haplotypes reconstructed based on Mendelian segregation rules were considered the gold standard to assess the performance of population-based methods in two scenarios. In the first one, only these 98 individuals were phased, while in the second one, all the 264 sequenced individuals were phased simultaneously, ignoring the pedigree relationships. We assessed phasing accuracy based on switch error counts (SEC) and rates (SER), lengths of correctly phased haplotypes and the probability that there is no phasing error between a pair of SNPs as a function of their distance. For most evaluated metrics or scenarios, the best software was either ShapeIT4.1 or Beagle5.2, both methods resulting in particularly high phasing accuracies. For instance, ShapeIT4.1 achieved a median SEC of 50 per individual and a mean haplotype block length of 24.1 Mb (scenario 2). These statistics are remarkable since the methods were evaluated with a map of 8,400,000 SNPs, and this corresponds to only one switch error every 40,000 phased informative markers. When more relatives were included in the data (scenario 2), FImpute3.0 reconstructed extremely long segments without errors. CONCLUSIONS: We report extremely high phasing accuracies in a typical livestock sample. ShapeIT4.1 and Beagle5.2 proved to be the most accurate, particularly for phasing long segments and in the first scenario. Nevertheless, most tools achieved high accuracy at short distances and would be suitable for applications requiring only local haplotypes.

The virulome of Streptomyces scabiei in response to cello-oligosaccharide elicitors.

The development of spots or lesions symptomatic of common scab on root and tuber crops is caused by few pathogenic Streptomyces with Streptomyces scabiei 87-22 as the model species. Thaxtomin phytotoxins are the primary virulence determinants, mainly acting by impairing cellulose synthesis, and their production in S. scabiei is in turn boosted by cello-oligosaccharides released from host plants. In this work we aimed to determine which molecules and which biosynthetic gene clusters (BGCs) of the specialized metabolism of S. scabiei 87-22 show a production and/or a transcriptional response to cello-oligosaccharides. Comparative metabolomic analyses revealed that molecules of the virulome of S. scabiei induced by cellobiose and cellotriose include (i) thaxtomin and concanamycin phytotoxins, (ii) desferrioxamines, scabichelin and turgichelin siderophores in order to acquire iron essential for housekeeping functions, (iii) ectoine for protection against osmotic shock once inside the host, and (iv) bottromycin and concanamycin antimicrobials possibly to prevent other microorganisms from colonizing the same niche. Importantly, both cello-oligosaccharides reduced the production of the spore germination inhibitors germicidins thereby giving the 'green light' to escape dormancy and trigger the onset of the pathogenic lifestyle. For most metabolites - either with induced or reduced production - cellotriose was revealed to be a slightly stronger elicitor compared to cellobiose, supporting an earlier hypothesis which suggested the trisaccharide was the real trigger for virulence released from the plant cell wall through the action of thaxtomins. Interestingly, except for thaxtomins, none of these BGCs' expression seems to be under direct control of the cellulose utilization repressor CebR suggesting the existence of a yet unknown mechanism for switching on the virulome. Finally, a transcriptomic analysis revealed nine additional cryptic BGCs that have their expression awakened by cello-oligosaccharides, suggesting that other and yet to be discovered metabolites could be part of the virulome of S. scabiei.

Repetitive saliva-based mass screening as a tool for controlling SARS-CoV-2 transmission in nursing homes.

Nursing home (NH) residents and staff have been severely affected by the COVID-19 pandemic. The aim of this study was to examine the use of weekly saliva RT-qPCR testing for SARS-CoV-2 detection among NH workers as a strategy to control disease transmission within NHs in Belgium. From 16 November to 27 December 2020, a voluntary and anonymous weekly screening was implemented in a cohort of 50,000 workers across 572 NHs in the Walloon region of Belgium to detect asymptomatic cases of SARS-CoV-2 via saliva RT-qPCR testing and using the Diagenode saliva sample collection device. Positive workers were isolated to avoid subsequent infections in residents and other staff. RT-qPCR testing was based on pooled saliva sampling techniques from three workers, followed by individual testing of each positive or inconclusive pool. The majority of NHs (85%) and 55% of their workers participated. Pooling did not affect sensitivity as it only induced a very decrease in sensitivity estimated as 0.33%. Significant decreases in the prevalence (34.4-13.4%) and incidence of NHs with either single (13.8-2%) or multiple positive workers (3.7-0%) were observed over time. In addition, deaths among NH residents and NH worker absences decreased significantly over time. Weekly saliva RT-qPCR testing for SARS-CoV-2 demonstrated large-scale feasibility and efficacy in disrupting the chain of transmission. Implementation of this testing strategy in NHs could also be extended to other settings with the aim to control viral transmission for maintaining essential activities.

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle.

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.

Single-cell RNA sequencing of human liver reveals hepatic stellate cell heterogeneity.

BACKGROUND & AIMS: The multiple vital functions of the human liver are performed by highly specialised parenchymal and non-parenchymal cells organised in complex collaborative sinusoidal units. Although crucial for homeostasis, the cellular make-up of the human liver remains to be fully elucidated. Here, single-cell RNA-sequencing was used to unravel the heterogeneity of human liver cells, in particular of hepatocytes (HEPs) and hepatic stellate cells (HSCs). METHOD: The transcriptome of ~25,000 freshly isolated human liver cells was profiled using droplet-based RNA-sequencing. Recently published data sets and RNA in situ hybridisation were integrated to validate and locate newly identified cell populations. RESULTS: In total, 22 cell populations were annotated that reflected the heterogeneity of human parenchymal and non-parenchymal liver cells. More than 20,000 HEPs were ordered along the portocentral axis to confirm known, and reveal previously undescribed, zonated liver functions. The existence of 2 subpopulations of human HSCs with unique gene expression signatures and distinct intralobular localisation was revealed (i.e. portal and central vein-concentrated GPC3 + HSCs and perisinusoidally located DBH + HSCs). In particular, these data suggest that, although both subpopulations collaborate in the production and organisation of extracellular matrix, GPC3 + HSCs specifically express genes involved in the metabolism of glycosaminoglycans, whereas DBH + HSCs display a gene signature that is reminiscent of antigen-presenting cells. CONCLUSIONS: This study highlights metabolic zonation as a key determinant of HEP transcriptomic heterogeneity and, for the first time, outlines the existence of heterogeneous HSC subpopulations in the human liver. These findings call for further research on the functional implications of liver cell heterogeneity in health and disease. LAY SUMMARY: This study resolves the cellular landscape of the human liver in an unbiased manner and at high resolution to provide new insights into human liver cell biology. The results highlight the physiological heterogeneity of human hepatic stellate cells.

Alzheimer's disease genetic risk and sleep phenotypes in healthy young men: association with more slow waves and daytime sleepiness.

STUDY OBJECTIVES: Sleep disturbances and genetic variants have been identified as risk factors for Alzheimer's disease (AD). Our goal was to assess whether genome-wide polygenic risk scores (PRS) for AD associate with sleep phenotypes in young adults, decades before typical AD symptom onset. METHODS: We computed whole-genome PRS for AD and extensively phenotyped sleep under different sleep conditions, including baseline sleep, recovery sleep following sleep deprivation, and extended sleep opportunity, in a carefully selected homogenous sample of 363 healthy young men (22.1 years ± 2.7) devoid of sleep and cognitive disorders. RESULTS: AD PRS was associated with more slow-wave energy, that is, the cumulated power in the 0.5-4 Hz EEG band, a marker of sleep need, during habitual sleep and following sleep loss, and potentially with larger slow-wave sleep rebound following sleep deprivation. Furthermore, higher AD PRS was correlated with higher habitual daytime sleepiness. CONCLUSIONS: These results imply that sleep features may be associated with AD liability in young adults, when current AD biomarkers are typically negative, and support the notion that quantifying sleep alterations may be useful in assessing the risk for developing AD.

Analyses of Avascular Mutants Reveal Unique Transcriptomic Signature of Non-conventional Endothelial Cells.

Endothelial cells appear to emerge from diverse progenitors. However, to which extent their developmental origin contributes to define their cellular and molecular characteristics remains largely unknown. Here, we report that a subset of endothelial cells that emerge from the tailbud possess unique molecular characteristics that set them apart from stereotypical lateral plate mesoderm (LPM)-derived endothelial cells. Lineage tracing shows that these tailbud-derived endothelial cells arise at mid-somitogenesis stages, and surprisingly do not require Npas4l or Etsrp function, indicating that they have distinct spatiotemporal origins and are regulated by distinct molecular mechanisms. Microarray and single cell RNA-seq analyses reveal that somitogenesis- and neurogenesis-associated transcripts are over-represented in these tailbud-derived endothelial cells, suggesting that they possess a unique transcriptomic signature. Taken together, our results further reveal the diversity of endothelial cells with respect to their developmental origin and molecular properties, and provide compelling evidence that the molecular characteristics of endothelial cells may reflect their distinct developmental history.

SNP-based quantitative deconvolution of biological mixtures: application to the detection of cows with subclinical mastitis by whole-genome sequencing of tank milk.

Biological products of importance in food (e.g., milk) and medical (e.g., donor blood-derived products) sciences often correspond to mixtures of samples contributed by multiple individuals. Identifying which individuals contributed to the mixture and in what proportions may be of interest in several circumstances. We herein present a method that allows to do this by shallow whole-genome sequencing of the DNA in mixed samples from hundreds of donors. We show the efficacy of the approach for the detection of cows with subclinical mastitis by analysis of farms' tank mixtures containing milk from as many as 500 cows.

Detection and identification of transgenic events by next generation sequencing combined with enrichment technologies.

Next generation sequencing (NGS) is a promising tool for analysing the quality and safety of food and feed products. The detection and identification of genetically modified organisms (GMOs) is complex, as the diversity of transgenic events and types of structural elements introduced in plants continue to increase. In this paper, we show how a strategy that combines enrichment technologies with NGS can be used to detect a large panel of structural elements and partially or completely reconstruct the new sequence inserted into the plant genome in a single analysis, even at low GMO percentages. The strategy of enriching sequences of interest makes the approach applicable even to mixed products, which was not possible before due to insufficient coverage of the different genomes present. This approach is also the first step towards a more complete characterisation of agrifood products in a single analysis.